The boundary between long-range and short-range critical behavior 
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We investigate phase transitions of two-dimensional Ising models with power-law interactions, 
using an efficient Monte Carlo algorithm. For slow decay, the transition is of the mean-field type; for 
fast decay, it belongs to the short-range Ising universality class. We focus on the intermediate range, 
where the critical exponents depend continuously on the power law. We find that the boundary 
with short-range critical behavior occurs for interactions depending on distance r as r~ 15 / 4 . This 
answers a long-standing controversy between mutually conflicting renormalization-group analyses. 
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PACS numbers: 05.70.Jk, 64.60.Ak, 64.60.Fr 

Many of the intermolecular forces that play a central 
role in large areas of chemistry, physics, and biology have 
a long-ranged nature. Well-known examples are electro- 
static interactions, polarization forces, and van der Waals 
forces. Remarkably, there are still considerable deficien- 
cies in our knowledge of the critical behavior induced by 
these interactions, which include the prominent case of 
Coulombic criticality (see Ref. jl| for a review). 

But also for the simpler case of purely attractive, long- 
range interactions, there exists a smothering controversy, 
which we aim to resolve in this work. The present un- 
derstanding of critical behavior in systems with such al- 
gebraically decaying interactions is largely based on the 
renormalization-group (RG) calculations for the O(n) 
model by Fisher, Ma, and Nickel ||. Their analysis 
revealed that different regimes can be identified for the 
universal critical properties, characterized by the decay 
power of the interactions. In view of the small number 
of global parameters determining the universality class, 
the location of the boundaries between these regimes is of 
considerable interest. It is, therefore, disturbing to note 
that there appears to be still no consensus on the theo- 
retical side regarding the precise location of the bound- 
ary between short-range and long-range critical behavior. 
Here, we address this issue by means of numerical calcu- 
lations, which allow us, with minimal prior assumptions, 
to decide between contradictory RG scenarios. 

Our approach specializes to the Ising model, n = 1 , in 
d dimensions, described by the reduced Hamiltonian 
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where the spins Sk — ±1 are labeled by the lattice site 
fc, the sum runs over all spin pairs, and the pair coupling 
depends on the distance = \ fi — fj | between the spins. 
According to the analysis of Fisher et al. , universality 
classes are parametrized by a, and the following three 
distinct regimes were identified: a) The classical regime. 
The upper critical dimension is given by d u — 2a, so 



that mean-field-type critical behavior occurs for a < d/2. 
(b) The intermediate regime d/2 < a < 2; here the crit- 
ical exponents are continuous functions of a. (c) The 
short-range regime. For a > 2 the universal properties 
are those of the model with short-range interactions, e.g., 
between nearest-neighbors only; thus one observes that 
for d = 3 van der Waals interactions (decaying as l/r e ) 
actually lie quite close to the boundary between regimes 
(b) and (c). 

Although the general outline of these results has been 
widely accepted, one part of this picture has become the 
scene of a debate that appears, up to now, to be un- 
resolved. This concerns the situation close to a = 2. 
In Ref. pi is has been conjectured that, throughout the 
intermediate regime (b), the correlation-function expo- 
nent r\ is exactly given by 2 — a. On the other hand, 
in the short-range regime (c), i] takes a constant (but in- 
dependent) value ?y sr > for all d < 4, resulting in a jump 
discontinuity in rj at the decay power a — 2. While this 
remarkable phenomenon is not forbidden by thermody- 
namic arguments (which only require r\ < 2 + a) , it has 
attracted considerable attention over the past decades 
and various efforts have been undertaken to reinvesti- 
gate the corresponding RG scenario (|, |J, H> Hi- I n ad- 
dition, it may be remarked that this scenario does not 
capture the one-dimensional case, where a phase tran- 
sition is rigorously known to be absent j7| for a > 1, 
rather than for a > 2 (see also Ref. || and references 
therein). The issue was first addressed by Sak ||, who 
pointed out that higher-order terms in the RG equations 
considered in Ref. [|| generate additional short-range in- 
teractions in the renormalization process, which affects 
the competition between the long-range and the short- 
range fixed points of the RG transformation. As a con- 
sequence, for d < 4 the boundary between the interme- 
diate and the short-range regime was found to shift from 
(T = 2toer = 2 — fj, where the e expansion for fj agrees 
to lowest orders with that of the short-range exponent 
7] SI . In a field-theoretic approach, Honkonen and Nal- 
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imov H proved, to all orders in perturbation theory, the 
stability of the short-range fixed point for a > 2 — r) SI 
and of its long-range counterpart for a < 2 — rfc, where 
r;i,. is the anomalous dimension of the field, evaluated at 
the long-range fixed point. Note that these authors also 
pointed out that the former result can be obtained from 
simple scaling arguments, but the latter cannot. Appeal- 
ing aspects of these findings are firstly the continuous and 
monotonic a dependence of the correlation-function ex- 
ponent (provided that r]\ T and 77 sr coincide at a — 2 — 7y sr ) , 
and secondly that the theory has now attained consis- 
tency with the exact results for the one-dimensional case. 

However, the analysis of || has also been the subject of 
criticism. Van Enter [H pointed out an apparent incon- 
sistency in Sak's results. Three-dimensional short-range 
O(n) models exhibit a nonzero spontaneous magnetiza- 
tion below the critical temperature. Van Enter || de- 
rived that, for n > 2, this broken symmetry is incon- 
sistent with the presence of a long-range perturbation 
with a < 2. This contradicts Sak's finding that such 
perturbations are irrelevant in a non-vanishing range 
2 — rjsr < a < 2. While this inconsistency emerges for 
n > 2 and antiferromagnetic perturbations only, it may 
be taken as a sign that the renormalization scheme |], |(| 
is incomplete, which could thus also affect the purely 
ferromagnetic n = 1 case. Another conflict with Sak's 
results arose from the analyses of Yamazaki Q . A qual- 
itative difference between his and Sak's result for the 
correlation-function exponent as a function of a is the 
absence of a kink at a = 2 — rj sr . However, we note a 
problem with the internal consistency of his results: The 
smoothly varying rj{a) was obtained by using a long- 
range (i.e., fc^-type, where k is the wave vector used 
in the Fourier representation) instead of a short-range 
(fc 2 -type) propagator in the RG calculations, despite the 
fact that, in the same calculation, it was found that the 
long-range term in the Landau- Ginzburg- Wilson Hamil- 
tonian is irrelevant, and only the short-range term sur- 
vives. A more serious objection was raised by Gusmao 
and Theumann || (but ignored in Ref. ||), who argued 
that the parameter 2 — er (i.e., 2a in the notation of 
Rcf. ||) is not a valid expansion parameter. They re- 
considered the problem using renormalizcd perturbation 
theory in s' — 2er — d, and found that the long-range 
expansion in e' is stable with respect to short-range per- 
turbations up to a = 2. This implies a restoration of 
the early result of Fisher et al. Q that the boundary be- 
tween the intermediate and the short-range regime lies at 
(7 = 2. While, at least for d = 3, the difference between 
this value and er = 2 — ?7 sr may be small, it indicates a 
fundamental dichotomy, which is also relevant for other 
systems with competing fixed points. 

Until recently, it appeared not to be feasible to shed 
some light on this unsatisfactory situation by means of 
Monte Carlo simulations, for the following reasons. Since 
the problems reside in the transition region between the 



intermediate and the short-range regime, one expects cor- 
rections to scaling to converge only slowly with increas- 
ing system size. This necessitates the simulation of large 
systems. Since each spin interacts with all other spins in 
the system, the amount of work per Metropolis-type spin 
visit then becomes prohibitive. Use of a cluster Monte 
Carlo algorithm JhJ that employs a number of opera- 
tions per spin flip that is independent of the system size, 
and moreover suppresses critical slowing down, has now 
enabled us to obtain highly accurate data for sufficiently 
large system sizes. 

We present numerical results for the critical exponent ij 
and the Binder cumulant JlT| as a function of er, that al- 
low us to distinguish between the different theories. The 
simulated systems are defined on L x L lattices with pe- 
riodic boundaries and sizes from L — 4 to L = 1 000. We 
chose to study two-dimensional systems not only to max- 
imize the attainable linear system size, but in particular 
because there the exponent rj sx = j has a much larger 
value than for d — 3 (r/ sr ~ 0.037); this maximizes both 
the size of the disputed region (2 — ry sr , 2) and the mag- 
nitude of the supposed jump in r)(o~). The length of the 
simulations was such, that (for the largest system sizes) 
typically a relative uncertainty of one part in a thousand 
was reached for the Binder cumulant. The precise form 
of the pair interaction was taken as 

m) = K Cl dx Ci dy (* 2 + ' (2) 

where f = (r x ,r y ) denotes the difference between the 
integer coordinates of the two interacting spins. It is im- 
portant to note that this differs from the interaction in 
Eq. (0) only in powers of r that decay faster than r — d— °" ; 
these terms are irrelevant and will only affect nonuniver- 
sal terms, such as the location of the critical point. The 
critical exponents and boundaries of the regimes (a)-(c) 
will not be altered. As mentioned, we have adopted pe- 
riodic boundary conditions, which are expected to speed 
up the approach to the thermodynamic limit. No artifi- 
cial cutoff was imposed on the interaction range: Inter- 
actions take place between all periodic images. 

We have carried out simulations for decay powers 
a = 1.2, 1.4, 1.6, 1.75, 1.85, 1.95, 2.0, 2.05, 2.25, 2.5, 
2.75, and 3.0. For each value of a, we have determined 
the critical coupling K c on the basis of a finite-size scal- 
ing analysis of the susceptibility-like quantity (to 2 ) and 
of the dimensionless amplitude ratio Q = (m 2 ) 2 /(m 4 ), 
where to is the magnetization density. Q is a variant of 
the fourth-order cumulant introduced by Binder p"l| . For 
Ising-like systems above the upper critical dimension, i.e., 
in the classical regime, it has been predicted ]l2j to take 
the universal critical value 87r 2 /r 4 (i) ~ 0.456947, the 
same value as in the mean-field model |l0| ; this has been 
confirmed numerically in Ref. p3| . In the short-range 
regime it takes a universal value close to Q = 0.856216 • • • 
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TABLE I: Results of detailed finite-size-scaling analyses for 
two-dimensional systems with power-law interactions, for sev- 
eral values of the decay parameter a. The results for the 
critical point K c , the fourth-order amplitude ratio Q and the 
correlation-function exponent 77 are followed by the estimated 
error in the last decimal places. The last column lists the 
predicted values for rj, see the text. 
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fu] , |Tj|. In the intermediate regime, finally, Q has been 
calculated by means of a singular expansion in e', up to 
second order in \fF ]l6], 0. From the magnetic suscep- 
tibility, which diverges as L^l v at criticality, r\ can be 
extracted using the scaling law ■) = (2 — //);/. An impor- 
tant issue in these finite-size analyses is the fact that the 
correction-to-scaling exponents are essentially unknown, 
and actually depend on the borderline value of a. In par- 
ticular, there will be logarithmically decaying finite-size 
corrections at the borderline itself. Considerable effort 
has been exercised to encompass all such uncertainties in 
the quoted margins for our estimates for K c , Q, and r\ 
(see Table |). 

While the detailed least-squares analysis depended on 
the value of cr, the following general procedure was ap- 
plied. An (effective) leading correction-to-scaling expo- 
nent yi was obtained from a fit of the amplitude ratio Q, 
in which also temperature-dependent terms were taken 
into account. Subsequently, 7/1/ was extracted from an 
analysis of (m 2 ) in which yi was kept fixed. This al- 
lowed us to include higher-order corrections. The anal- 
ysis was then repeated with t/, varied within its uncer- 
tainty margins, in order to determine the uncertainty in 
7/1/. Finally, the entire analysis was repeated for differ- 
ent minimum system sizes. For a = 1.75, a consistent 
analysis could only be obtained by allowing for logarith- 
mic corrections to scaling, as would be expected if the 
crossover from long-range to short-range criticality oc- 
curs at a — 2 — rj sx . 

The critical coupling exhibits a very smooth but also 
remarkable, almost parabola-shaped dependence on <r, 
with a maximum around a — 2.17. This nonmonotonic 
variation may be due to our use of the integrated cou- 
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FIG. 1: The amplitude ratio Q as a function of the decay 
parameter a. The behavior is consistent with a crossover from 
long-range to short-range critical behavior around a = 1.75, 
as predicted in Refs. fel 6[. 



pling (g) which tends to enhance the close-neighbor cou- 
plings, in particular for large a. For a > 2, both Q and 
77 are fully consistent with the known values for the two- 
dimensional Ising universality class. On the other hand, 
for a < 1.60, 77 is in excellent agreement with the RG pre- 
diction 2 — o ||] ; Q exhibits a remarkably linear deviation 
from its mean-field value ]l6]] . However, in the disputed 
region 1.75 < a < 2.0, where we invested most of our 
computing efforts, the uncertainties are still considerable; 
this is by itself already indicative of the occurrence of a 
crossover between two fixed points. 

In the absence of accurate predictions for the ampli- 
tude ratio Q (depicted in Fig. |l]), this quantity does not 
provide much further evidence to distinguish between the 
different theoretical scenarios. However, there is no evi- 
dence for a jump in Q at a = 2; this fails to support the 
Gusmao and Theumann H analysis. For a = 1.85 and 
(T = 1.95, we find that Q is compatible with the 2D Ising 
value, and the deviation of about 1.5 standard errors at 
(j = 1.75 is still within acceptable limits. 

On the other hand, the exponent 77 offers clear evidence 
in order to distinguish between the different scenarios, cf. 
Fig. ||. While also this quantity suffers from relatively 
large uncertainties in the crossover region, these uncer- 
tainties are small enough to exclude the scenario of Ref. 
H which implies that a transition to short-range critical 
behavior occurs at a — 2 and that rj = 2— a for a < 2. In- 
deed, instead of the 2D Ising value r\ = j , this exponent 
should then take the values 0.10 at a — 1.85 and 0.05 
at a = 1.95. However, at these two values of cr, our re- 
sults for rj lie within one standard deviation from the 2D 
Ising value, and three respectively almost five standard 
deviations above 2 — cr. 

In summary, we have demonstrated that both the 
correlation-function exponent 77 and the fourth-order am- 
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the limit a — > 2, as is done in Ref. M. This would also 
explain why the long-range (e') expansion does not pro- 
vide any indication of a change to short-range behavior 
at a — 2 — rj ST . 

We are indebted to Professor J. Sak for providing his 
valuable comments on the status of the theoretical analy- 
ses of the border between the regions of the intermediate- 
and the short-range interactions. This research is sup- 
ported in part by the Dutch FOM foundation ("Sticht- 
ing voor Fundamenteel Onderzoek der Materie") which 
is financially supported by the NWO ( "Nederlandse Or- 
ganisatie voor Wetenschappelijk Onderzoek"). 
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FIG. 2: The correlation-function exponent r\ as a function 
of the decay parameter a. For a < 1.75 the data are in 
excellent agreement with the conjecture r\ — 2 — a Ej . whereas 
for <F~ 
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plitude ratio Q take their (universal) short-range Ising 
values for a > 2 — J] SI . For a > 2 this could be shown 
with high numerical accuracy, but also for 2 — ?7 sr < a < 2 
the results are precise enough to convincingly exclude 
a crossover from long-range to short-range critical be- 
havior at a = 2. Instead, they support a crossover at 
(7 = 2 — r/ sr , in agreement with Refs. [|[ This pro- 
vides a clear answer to the dilemma posed by conflicting 
renormalization-group scenarios. 

Furthermore, our findings for r\ are in excellent agree- 
ment with r] — 2 — a in the intermediate range d/2 < a < 
2 — i] sr , supporting the conjecture that all higher-order 
contributions vanish in the e' expansion for r\ || and in 
contrast with the scenario proposed in Ref. Q . 

The amplitude ratio Q depends approximately linearly 
on <7 for d/2 < a < 2 — r\ sr . Remarkably, the most promi- 
nent deviations from linearity occur near a = 2 — i] sr , 
while the e' expansion predicts a square-root-like singu- 
larity at the opposite end of the intermediate range [jl6| . 

Finally, we quote a plausible explanation for the break- 
down of the description of Ref. Q in the range 2 — r/ sr < 
a < 2 jUI]. If the validity of the e'-expansion is not uni- 
form in a, the region of convergence may shrink to zero 
when a approaches 2. Then one cannot fix e' and take 
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